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Abstract 

1 + 

We present results for the first positive parity excited state of the nucleon, namely, the Roper resonance (N^ =1440 
MeV) from a variational analysis technique. The analysis is performed for pion masses as low as 224 MeV in 
quenched QCD with the FLIC fermion action. A wide variety of smeared-smeared correlation functions are used to 
construct correlation matrices. This is done in order to find a suitable basis of operators for the variational analysis 
such that eigenstates of the QCD Hamiltonian may be isolated. A lower lying Roper state is observed that approaches 
the physical Roper state. To the best of our knowledge, the first time this state has been identified at light quark masses 
using a variational approach. 
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One of the long-standing puzzles in hadron spec- 
troscopy has been the low mass of the first positive par- 
ity, - J , excitation of the nucleon, known as the 
Roper resonance A^*(1440 MeV). In constituent or va- 
lence quark models with harmonic oscillator potentials, 
the lowest-lying odd parity state naturally occurs be- 
low the N = J state (with principal quantum number 
N - 2) ili 01 whereas, in nature the Roper resonance 
is almost 100 MeV below the N = i"(1535 MeV) state. 
Similar difficulties in the level orderings appear for the 
J'' = |^A*(1600) and i^2*(1690) resonances, which 
have led to the speculation that the Roper resonance 
may be more appropriately viewed as a hybrid baryon 
state with explicitly excited gluon field configurations 
islS]^ or as a breathing mode of the ground state [5] or 
states which can be described in terms of meson-baryon 
dynamics alone Jil]. 

The first detailed analysis of the positive parity exci- 
tation of the nucleon was performed in Ref. 101 using 
Wilson fermions and an operator product expansion 
spectral ansatz. Since then several attempts have been 
made to address these issues in the lattice framework 
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been identified using Bayesian techniques. 



II [3 [11 [H, H, but in many 



cases no potential identification of the Roper state has 



been made 
analysis of 




ll 1II12I1 . Recently however, in the 
llSll . a low-lying Roper state has 
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Here, we use a 'variational method' 2^ 21], 

which is based on a correlation matrix analysis and has 
been used quite extensively in Refs. ll ll llo. 2l[ 22, 

M [35, 



23 

36l l37[l . Though the ground state mass of the nucleon 
has been described successfully, an unambiguous deter- 
mination of the Roper state with this method has not 
been achieved in the past, though significant^ amounts 
of research have been carried out in Ref. |23[l . by the 
CSSM Lattice Collaboration HH [H iSTthe BGR 
Mr 27, 28, 29, 33] Collaboration and in Refs. HEl. 
In this Letter, we present evidence of a low-lying Roper 
state for the first time using a variational analysis. The 
observed state displays chiral curvature and approaches 
the physical mass of the Roper state. The standard nu- 
cleon interpolating field xi is considered in this analy- 
sis. Various sweeps of Gaussian smearing |38] are used 
to construct a smeared-smeared correlation function ba- 
sis from which we obtain the correlation matrices. 
The two point correlation function matrix for p - can 

August 26. 2009 



be written as 



z 



(2) 



where, Dirac indices are implicit. Here, Aj and Aj are 
the coupUngs of interpolators and^y at the sink and 
source, respectively, a enumerates the energy eigen- 
states with mass OTq,. 

Since the only t dependence comes from the exponen- 
tial term, one can seek a linear superposition of inter- 
polators, X iu", such that (more detail can be found in 
Refs. LiLi6]), 



Gij{t + M) u'] ^e-"'-''"Gij{t)u"j, 



(3) 



for sufficiently large t and t + At, see Refs. [22] and 
lfl6s1 . Multiplying the above equation by [G,)(f)]"' from 
the left leads to an eigenvalue equation. 



[iGit)r^G(t + At)]iju"j = c" u[ 



(4) 



where c" = e^"'"'^' is the eigenvalue. Similar to Eq.(|4]i, 
one can also solve the left eigenvalue equation to re- 
cover the v" eigenvector. 



vf [G(f + Af)(G(f)r'],v = c%°'. 



(5) 



The vectors u" and v" diagonalize the correlation matrix 
at time t and t + At making the projected correlation 
matrix. 



vfG,/f)«^ 5"". (6) 
The parity projected, eigenstate projected correlator. 



vfGf/f)«J ^ Gl, 



(7) 



is then analyzed using standard techniques to obtain 
masses of different states. 

Our analysis is exploratory, seeking to develop tech- 
niques to access the Roper state in lattice gauge theory. 
Our lattice ensemble consists of 200 quenched configu- 
rations with a lattice volume of 16^ x 32. Gauge field 
configurations are generated by using the DBW2 gauge 
action ll39i l40ll and an C?(fl)-improved FLIC fermion 
action 14111 is used to generate quark propagators. This 
action has excellent scaling properties and provides near 
continuum results at finite lattice spacing |42]. The lat- 
tice spacing is a = 0.1273 fm, as determined by the 
static quark potential, with the scale set using the Som- 
mer scale, ro = 0.49 fm 143}]. In the irrelevant op- 
erators of the fermion action we apply four sweeps of 



stout-link smearing to the gauge links to reduce the cou- 
pling with the high frequency modes of the theory ll44ll 
providing 0{a) improvement We use the same 

method as in Ref. [16, 45J to determine fixed bound- 
ary efifects, and the eff'ects are significant only after 
time slice 25 in the present analysis. Various sweeps 
(1, 3, 7, 12, 16, 26, 35, 48 sweeps corresponding to 
rms radii, in lattice units, of 0.6897, 1.0459, 1.5831, 
2.0639, 2.3792, 3.0284, 3.5237, 4.1868) of gauge in- 
variant Gaussian smearing [38] are applied symmetri- 
cally at the source (at t - A) and at the sink. This is to 
ensure a variety of overlaps of the interpolators with the 
lower-lying states. The analysis is performed on ten dif- 
ferent quark masses corresponding to pion masses m„ - 
{0.797, 0.729, 0.641, 0.541, 0.430, 0.380, 0.327, 0.295, 
0.249,0.224} GeV. Error analysis is performed using 
a second-order single elimination jackknife method, 
where the ;tf^/dof is obtained via a covariance matrix 
analysis method. Our fitting method is discussed exten- 
sively in Ref. [[Ifl. 

The nucleon interpolator we consider is the local scalar- 
di quar k interpolator having a non-relativistic reduction 



Xi(x) = e"''\u"'{x)Cysd\x))u%x). 



(8) 



We consider several 4x4 matrices. Each matrix is 
constructed with diff'erent sets of correlation functions, 
each set element corresponding to a different numbers 
of sweeps of gauge invariant Gaussian smearing at the 
source and sink of the X\X\ correlators. This provides 
a large basis of operators with varieties of overlap 
among energy states. We consider seven combina- 
tions {1=(1,7,16,35), 2=(3,7,16,35), 3=(1, 12,26,48), 
4^(3,12,26,35), 5=(3,12,26,48), 6^(12,16,26,35), 
7=(7,I6,35,48)) of 4 x 4 matrices. In Ref. [16] it 
was shown that one cannot isolate a low-lying excited 
eigenstate using a single fixed-size source smearing. 
The superposition of states manifested itself as a 
smearing dependence of the effective mass. In this 
Letter we exploit this sensitivity to isolate the energy 
eigenstates. 

In Figdl we show the mass from the projected cor- 
relation functions and from the eigenvalues (as shown 
in Ref. (Il6l|) for the fourth combination (3,12,26,35) of 
4x4 matrices. We note that similar results in mass from 
the projected correlation functions and the eigenvalues 
are observed in this analysis as in Ref. m. Though 
the mass of the first excited state from projected corre- 
lation functions show little change with variational pa- 
rameters, see FiglT] the second and third excited states 
start a little below which indicates t and t+At are not suf- 
ficiently large. With larger Euclidean times fewer states 
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Figure 1: (Color online). Masses of the nucleon, A^^ - states, from projected correlation functions as shown in Eq.JT) (left figure) and from 
eigenvalues (right figure), for the pion mass of 797 MeV, and for the 4* combination (3,12,26,35) of 4 x4 matrices. Each set of ground and excited 
states masses correspond to the diagonalization of the correlation matrix for each set of variational parameters t = tstmt (shown in major tick marks) 
and a; (shown in minor tick marks). In the legend "g.s" stands for the ground state, whereas, "e.s" is for excited state. Larger values of tstan and At 
did not provide a stable eigenvalue analysis. 
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Figure 2: (Color online). Masses of the nucleon, A'l - states, from projected correlation functions as shown in Eq.0 for the pion mass of 797 
MeV. Numbers in the horizontal scale correspond to each combination of smeared 4x4 correlation matrices. For instance, 1 and 2 correspond to 
the combinations of (1,7,16,35) and (3,7,16,35) respectively and so on, as discussed in the text following Eq.jS). Masses are extracted according to 
the selection criteria described in Ref. Ilql. 
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will contribute significantly to the correlators. The ro- 
bust aspect of fitting projected correlators is manifest 
in Figd] and reflects the stability of the eigenvectors 
against changes in t and t + At. In contrast, the mass 
from the eigenvalue analysis shows significant depen- 
dence on the variational parameters. The same method 
as described in Ref. [16] is applied in this Letter to ex- 
tract the mass from the projected correlation functions. 

In Fig|2] masses extracted from all the combinations 
of 4 X 4 matrices (from to 7*) are shown for the 
pion mass of 797 MeV. Some dependence of the ex- 
cited states on smearing count is also observed here as 
in Ref. [16] for a few of the interpolator basis smearing 
sets. However the ground and first excited states are ro- 
bust against changes in the interpolator basis, providing 
evidence that an energy eigenstate has been isolated. It 
should be noted that the highest excited state (the third 
excited state) is influenced more by the level of smear- 
ing than the lower excited states. This is to be expected 
as this state must accommodate all remaining spectral 
strength. 

The P' combination in Fig|2]provides heavier excited 
states as this basis begins with a low number of smear- 
ing sweeps (a sweep count of 1) and also contains an- 
other low smearing set of 7 sweeps. The second and 
third excited states, and more importantly, the first ex- 
cited state sits a bit high in comparison with the other 
bases. Hence, extracting masses with this basis is not as 
reliable as other sets. The 2°'* combination also contains 
elements with a small smearing sweep count (3 and 7), 
hence this basis also provides heavier excited states and 
shows some systematic drift in the second excited state. 
However, this basis has reduced contamination from the 
excited states when compared with the first basis. The 
3"* combination also starts at the low smearing count, 
so the mass from this basis for the third excited state is 
a little high. 

We can observe at this point that including basis el- 
ements with a low smearing count will increase the 
masses of excited states (for instance, consecutive low 
numbers of smearing sweeps 1,7 and 3,7, respectively). 
This is because the correlation functions with these low 
sweep counts have a large overlap with several heavier 
excited states in their sub-leading exponential. We also 
observe that the inclusion of basis elements with a high 
level of smearing (for instance, a sweep count of 48) 
results in larger statistical errors in the analysis. 

The 6* combination starts from a moderate smear- 
ing sweep count (12) and also contains elements with 
consecutive smearing sweep counts which provides less 
diversity in the basis. Similarly, the 7* combination 
contains consecutive large smearing sweep counts of 35 



and 48 and operators for these levels of smearing are 
very similar challenging the isolation of single energy 
eigenstates. 

The 4* and 5* combinations are well spread over 
the given range of smearing sweeps and start from a 
sweep count of 3. They don't include successive lower 
smearing sweep counts. The 5* combination contains 
the basis element with a sweep count of 48 but has 
only slightly larger statistical errors than the 4* basis 
choice. Both these bases provide diversity. It is ob- 
served that both combinations provide consistent results 
for the states. Nonetheless, it should be noted that the 
^rd gth combinations all provide very consistent re- 
sults for the lower three (ground, first and second) en- 
ergy states, shown in Fig|2] While the 2"'^ and the 
3'^'' excited states display some evidence of eigenstate 
mixing, the ground and the excited states are ro- 
bust in this analysis as they agree within one standard 
deviation for all the combinations of 4 x 4 matrices. 
Hence, an analysis is performed to calculate the sys- 
tematic errors associated with the choice of basis over 
the preferred four combinations (from 3"* to 6*) with 

CTb - -^Jjjpi llfj^iiMi - M)2, where, A^^, is the number 
of bases, in this case is equal to 4. It should be noted 
that for the ground and first excited states the systematic 
errors associated with choice of basis are very small in 
comparison with their statistical errors (see Table[T)- 

In Figl3] the masses from projected correlation func- 
tions are shown for all pion masses averaged over four 
combinations of correlation matrices (from 3'^'' to 6*). 
Masses are averaged over these four bases and errors 
(average statistical errors over these four bases and sys- 
tematic errors associated with basis choices) are com- 
bined in quadrature, cr = ^cr^ + cr^. Eigenstate mixing 
may be affecting the results for the second and the third 
excited states. The ground and the I'*' excited states sys- 
tematically approach the physical mass of the Nucleon 
and Roper state. This is the first time evidence of a low- 
lying Roper state has appeared from variational analysis 
near the chiral regime. 

In FiglH the ground and the first excited states of 
Fig|3] are shown in larger scale for clarity. The non- 
interacting two particle P-wave N + n is shown by the 
dashed line. It is interesting to note that the observed 
lattice Roper state sits lower than the P-wave N + n, 
indicative of attractive ;rN interactions producing a res- 
onance at physical quark masses. The very consistent 
Euclidean-time fit window for all the bases and for all 
the quark masses, noticeably the lower time of the fit 
as shown in Table |2] assist to identify the robust nature 
of the extracted Roper state which comes at the same 
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Figure 3: (Color online). Mass of the nucleon, N'^ - states, for the ground and the excited states. The eiTors shown in the figure are a combination 
of average statistical errors and systematic eiTors due to basis choices over four bases (from 3"^ to 6"" of Fig[2j. Errors are combined in quadrature. 
The black filled symbols are the experimental values of the ground and the Roper states of the nucleon |47]. The rightmost point corresponds to 
the pion mass of 797 MeV, then for 739, 641, 541, 430, 380, 327, 295, 249 and 224 MeV (leftmost point). 



Table 1: Mass of the nucleon, A^J - states, are averaged over the four bases (from 3"* to 6*) and the errors for the nucleon shown here are a 
combination of average statistical errors and systematic en'ors for choice of basis over the four bases. The 2"'' and 3"* states may be superposition 
of energy eigenstates as discussed in the text. 
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0.5141(19) 


1.0399(65) 


1.457(41) 


1.827(62) 


2.464(78) 


0.4705(20) 


0.9905(71) 
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0.9280(79) 
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0.783(11) 


1.239(51) 


1.73(12) 
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1.07(11) 


1.885(74) 


2.37(26) 


0.1448(44) 


0.661(15) 


1.05(13) 


1.885(73) 


2.53(16) 
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Table 2: Mass of the Roper state, , for four bases {(1,12,26,48), (3,12,26,35), (3,12,26,48) and (12,16,26,35)1 of 4 X 4 matrices, fi and t2 
represent the lower and upper times respectively of fit windows in the projected effective mass. Inclusion of ?i = 6 produces x^/'iof 3> I. 
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Figure 4: (Color online). The ground and the Roper states and the non-interacting P-wave N + n are illustrated. The black filled symbols are the 
experimental values of the ground and the Roper states obtained from Ref. BtIi . 
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Euclidean-time from the heaviest to the hghtest quarks 
masses. 

There are quenched artifacts that will have significant 
influence on our results should we progress to lighter 
quark masses. Not only are nN couplings different in 
quenched QCD, but there are also contributions from 
unphysical, degenerate P-wave rj'N two-particle states 
lll4ll . Future calculations approaching the decay thresh- 
old should be done in full dynamic al-fermion QCD. 

In conclusion, through the use of a variety of 
smeared-smeared correlation functions in constructing 
the correlation matrix, the first positive parity excited 
state of the nucleon A^^ , the Roper state, has been ob- 
served for the first time using the variational analysis. 
While the 3x3 correlation matrix analysis of standard 
interpolators is insufficient to isolate these excited en- 
ergy eigenstates of QCD [|l6t] . using the new correla- 
tion matrix construction with smeared-smeared correla- 
tors enables us to extract the otherwise elusive Roper 
state. 
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